Phase Coexistence of Complex Fluids in Shear Flow 
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We present some results of recent calculations of rigid rod-like particles in 
shear flow, based on the Doi model. This is an ideal model system for exhibiting 
£N| | the generic behavior of shear-thinning fluids (polymer solutions, wormlike micelles, 

surfactant solutions, liquid crystals) in shear flow. We present calculations of phase 
coexistence under shear among weakly-aligned (paranematic) and strongly-aligned 
phases, including alignment in the shear plane and in the vorticity direction (log- 
rolling). Phase coexistence is possible, in principle, under conditions of both com- 
mon shear stress and common strain rate, corresponding to different orientations of 
' the interface between phases. We discuss arguments for resolving this degeneracy. 

Calculation of phase coexistence relies on the presence of inhomogeneous terms in 
q . the dynamical equations of motion, which select the appropriate pair of coexisting 

states. We cast this condition in terms of an equivalent dynamical system, and ex- 
plore some aspects of how this differs from equilibrium phase coexistence. 
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<N ! L INTRODUCTION 

<N . 

Shear flow induces phase transitions and dynamic instabilities in many complex fluids, in- 
cluding wormlike micelles [|T]-0], liquid crystals, [|5|-|Tl~|]; and lamellar surfactant systems which 
can "roll" into multilamellar vesicles ("onions") Jl~2|-|T5|]. These instabilities typically manifest 
themselves in non-monotonic constitutive curves such as those in Fig. [l] [JT6|— pT]] , and in several 
systems, including wormlike micelles and lamellar surfactant solutions, are accompanied by ob- 
servable coexistence of two macroscopic "phases" of material. 

If a mean strain rate forces the system to lie on an unstable part of the flow curve (with negative 
slope), the system can phase separate into regions (bands) with high and low strain rates and 
still maintain the applied strain rate [J7UrMjJE2n. Fig. [T| shows that phase separation can occur at 
either common stress or common strain rate, depending on the geometry of the bands [|23|]: bands 
stacked along the axis of a Couette cell have the same strain rate and different shear stresses, 
while radial phase separation imposes a uniform shear stress and different strain rates. The shear- 
thinning wormlike micelle system phase separates radially into common-stress bands, while shear- 
thickening systems have been observed to separate into bands with both the common strain rate 



(worms) Q24| , |25| ] or common stress (worms and onions) Q26| , [27[ ] geometries, although the evidence 
for true steady state phase separation at common stress is not yet firm. 

Other systems with flow-induced "phase transitions" include colloidal suspensions of plate- 
like particles (which shear thicken and sometimes crystallize) [ |Z%| ] or nearly monodisperse spheres 
], as well as a variety of surfactant-like solutions of diblock copolymers in selective solvents 



[ftO[]. With so many increasingly detailed and careful experiments on so many systems, it would be 
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nice to have a consistent framework for non-equilibrium transitions. Unfortunately, most systems 
are sufficiently complicated that none of the observed transitions can be completely described, 
even qualitatively, by a credible microscopic model. For example, in certain limits the class of 
shear-thinning wormlike micelle systems which shear band has a mature theory for the linear 
rheology, which may be extended using successful ideas from entangled polymer dynamics to 
predict an instability [jT8|]. Unfortunately, a complete description of phase coexistence also requires 
knowledge of the shear-induced state, as well as details of the concentration dependence and, as 
we shall see, the inhomogeneous contribution to the dynamical equations. We are, at present, far 
from having all of these ingredients, and in many cases we do not have a clear understanding of 
even the structure of the high shear rate state, much less its dynamics. 




strain rate y 



FIG. 1. Stress-strain-rate curves for the Doi model (from Fig. |2|) for two concentrations. The dashed 
lines are unstable steady states. Straight lines indicate possible coexistence (at the same concentration) with 
either common stress (horizontal lines) or strain rate (vertical line) in the two phases. Inset: geometries for 
common stress (left) or common strain rate (right) coexistence. 

Recently, we have studied a well-known model, the Doi model for rigid rod suspensions in 
shear flow pl"|], which, while admittedly the product of many approximations, provides physically 
well-founded dynamics for both quiescent and shear-induced states. Although the shear rates nec- 
essary for inducing a transition are, in practice, typically quite high unless very long rods are used, 
and physical systems are often susceptible to various dynamical instabilities, this model system is 
quite helpful for building intuition about how to calculate non-equilibrium "phase diagrams" and 
how their resulting topologies resemble and differ from their equilibrium counterparts. 

A vexing question for non-equlibrium calculations is how to replace the free energy minimiza- 
tion familiar from equilibrium thermodynamics, to determine the analog of a first order phase 
transition. In the context of Fig. [j], one needs to determine the selected stress for an imposed strain 
rate (for phase separation at a common stress). It has emerged that an unambiguous resolution of 
this problem is to include explicit non-local terms in the dynamical equations and explicitly con- 
struct the coexisting state [Q[7|,^3|,[S3],[nj. This reduces to the equilibrium construction in the case 
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of zero shear, and can be shown to yield a single (barring accidental degeneracies) stress (given all 
other imposed conditions) at which coexistence occurs [|35|]. 

Below we summarize some results of our calculations on the Doi model for rigid rod suspen- 
sions in shear flow JT0|] ; the details will be published elsewhere [|3~6|]. This system is surprisingly 
rich, given its apparent simplicity and fairly obvious coupling of internal order to flow. Then we 
discuss some aspects of the interface construction for determining coexistence, and how it com- 
pares to its equilibrium counterpart. 



H. THE DOI MODEL 

The modified Doi model [ |3T|JTC| ] describes the dynamics of the rod-like particles suspension. 
The orientational degrees of freedom are parametrized by the conventional liquid crystalline order 
parameter tensor 

Qap{*) = {v a vp ~ §<W), (2.1) 

where (■) denotes an average around the point r of the second moment of the rod orientations v. 
For rigid rods the phase diagram, and in fact the dynamics, can be more conveniently represented 
by the exclude volume parameter u, defined by 

u = (pLa, (2.2) 

where <p is the rod volume fraction, L is the rod aspect ratio and a is an Oil) prefactor pl|]. 
Beginning from the Smoluchowski equation for a solution of rigid rods, and including a Maier- 
Sauper-like orientational interaction parameter, Doi was able to derive approximate coupled equa- 
tions of motion for the dynamics of Q, and the fluid velocity v(r), including the liquid crystalline 
contribution to the fluid stress tensor. 

The essential physics is that flow tends to align the rodlike molecules, typically roughly paral- 
lel to the flow direction, and hence stabilizes a nematic, or aligned, state. To study other complex 
fluids we would have a structural variable analogous to Q; e.g. in the wormlike micelle system 
we might need, in addition to the orientation tensor, the dynamics of the mean micellar size. We 
have augmented the Doi model by allowing for concentration diffusion driven by chemical poten- 
tial gradients, included the dynamical response to inhomogeneities in liquid crystalline order and 
concentration, and included the translational entropy of mixing which gives the system a biphasic 
coexistence regime in the absence of flow. For a given stress or strain rate we determine phase co- 
existence by explicitly constructing a stable coexisting steady state, which require inhomogeneous 
terms in the equations of motion (arising here from free energy terms which penalize inhomo- 
geneities in Q and 0). This procedure and the model have been documented elsewhere [f7|JT0|], and 
the interface construction will be discussed in more detail in Sec. 

To calculate the phase diagram we solve for the steady state homogeneous solutions to the 
coupled dynamical equations for {(ft, Q, V}. This yields a set of solutions which are then candi- 
dates for phase coexistence. Coexistence is possible with either common stress or common strain 
rate in the coexisting phases, depending on geometry, and must be examined for all pairs of stable 
homogeneous states. We expect coexisting states to have different values for (ft, Q, and either the 
stress or strain rate. As mentioned above, we determine coexistence by finding the locus of control 
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parameters for which a stable interfacial solution between two homogeneous states exists. We 
parametrize the shear stress and strain rate as 
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(2.3) 
(2.4) 



where D ro is the rotational diffusion coefficient, and v\ and u 2 are 0(1) geometric constants. 

The Doi model has three stable steady states in shear flow [ J57[ ] : A weakly-ordered paranematic 
state I, with the major axis of the order parameter in the shear plane; & flow-aligning state N, with 
a larger order parameter and major axis in the shear plane; and a log-rolling state L, with major 
axis in the vorticity direction. Fig. [Z] shows homogeneous constitutive relations for the I and 
N states. As can be seen, the N and L states are, successively less viscous than the I phase at 
the same concentration, with a viscosity which decreases slightly with increasing concentration 
(reflecting the greater order and hence lower viscosity of more concentrated phases), in contrast 
to the less-ordered I phases, whose viscosity increases with concentration, as is usual for colloidal 
suspensions. 
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FIG. 2. Homogeneous stress a xy vs. strain rate 7 curves for various excluded volumes u (u is propor- 
tional to (p) and L = 5. Shown are the I and N branches. The L branch is only stabilized for the high 
concentration (u = 2.7). 



A. Common Stress Phase Separation 

We first discuss the phase diagrams for common stress coexistence, in which the phase sepa- 
ration is radial in a cylindrical Couette flow. For common stress coexistence of two phases / and 
//, the fraction £ in phase I is determined by the lever rule, 

4> = C<t>! + (1 - C)0i/ (2.5a) 
7" = or, + (1 - C)T», (2-5b) 

where and 7 are mean values. Fig. |3| shows the phase diagram calculated for l-N coexistence 
for L = 5. The tie lines denoting pairs of coexisting phases are horizontal in the a xy — u plane, 
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and have positive slopes in the (7 — u) plane because the more concentrated nematic phase flows 
faster at a given stress. For weak stresses the equilibrium system is slightly perturbed and the tie 
lines are almost horizontal, while at high stresses the tie lines become steeper as the composition 
difference between the phases decreases and vanishes at a critical point. 




2.55 2.60 2.65 2.70 ' 2.55 2.60 2.65 



u u 
FIG. 3. Phase diagram in the (a xy — u) (a) and (7— u) (b) planes, along with limits of stability of I and 
N phases. Tie lines (horizontal in (a), sloped in (b) connect coexisting phases. 




mean strain rate 7 * 10 4 7 * 10 4 7 * 10 4 

FIG. 4. Stress a xy vs. mean strain rate 7 for common stress coexistence. Solid lines are constitutive 
curves for the I and N branches; dotted lines and o denote the N branch with which each I branch coexists at 
the low strain rate coexistence boundary. • denotes banded stresses, whose plateaus do not satisfy an equal 
area construction with the homogeneous constitutive curve. 

Mean Constitutive Relations — From the information in Fig. [3|b we can calculated the mean con- 
stitutive relation that would be measured in an experiment on a system with a given prescribed 
mean concentration. Upon applying stress at a given concentration, the system traces a vertical 
path through Fig. |] until the two-phase region is reached, during which a xy {^) varies smoothly. 
At this stress a tiny band of N phase develops, with composition and strain rate determined by 
the lever rule; and o- xy (y) is non-analytic (Fig. f|), exhibiting a change in slope. As the stress and, 
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hence the mean strain rate, increases further the system visits successive tie lines in the 7— u plane, 
each with a higher stress and mean strain rate and different coexisting concentrations. For close 
to the equilibrium l-N transition (Fig. |]c) the tie lines in the (7 — u) plane are fairly flat and the 
stress a xy changes significantly through the two-phase region. More dilute systems (Fig. f|a,b) 
have steeper tie lines and straighter and flatter 'plateaus' in a xy (j). 

Controlled strain rate experiments should follow the homogenous flow curves, except for the 
coexistence regime. In this case, analogy with equilibrium systems suggests that the system should 
eventually nucleate into a phase-separated banded state, with a corresponding stress change. Ex- 
periments on wormlike micelles display this kind of behavior upon increasing the strain rate above 
that of the phase boundary. If the mean strain rate is on an unstable part of the flow curve, we ex- 
pect a 'spinodal' (or mechanical) instability. There is a small region (inside the loop in Fig. |3]a) 
where the system is unstable when brought, at controlled strain rate, into this region from either 
the I or N states. This corresponds to constitutive curves with the shape of curve b in Fig. [Z| 

For controlled stress experiments, for stresses larger than the minimum coexistence stress and 
less than the local maximum we expect the system to follow the homogenous flow curve until a 
nucleation event occurs (Fig. |]b). Then, the strain rate should increase to either that of proper 
banded or single N phase strain rate, depending on the magnitude of the stress. For stresses larger 
than the I limit of stability we expect, again by analogy with equilibrium, a spinodal-type insta- 
bility. This simple picture is not quite corroborated in wormlike micelles: Grand et al. reported 
that a stress within a narrow range above the coexistence stress could be applied, and the system 
remained on the "metastable" branch indefinitely. 

Analogy with equilibrium systems suggests similar behavior upon reducing the stress or strain 
rate from the high-shear branch. Careful experiments (on any system) are needed to test the idea. 
For example, it is interesting to examine that, upon reducing the strain rate below the upper strain 
rate for the onset of shear banding and above the limit of stability of the high strain rate branch, 
whether the stress would spontaneously increase into the banded state. 

8 
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FIG. 5. Composite phase diagrams for l-L and l-N coexistence at common stress. This represents two 
overlayed phase diagrams, and not a single phase diagram. 

Log Rolling Phase — Fig. ^] shows the phase diagram for paranematic-log rolling l-L coexistence. 
For non-zero stress the biphasic region shifts to higher concentrations, since the stability limit of 
the L phase shifts to higher concentrations. Since the I and L phases have major axes of alignment 
in orthogonal directions, there is not a critical point; rather, the biphasic region ends when the I 
phase becomes unstable to the N phase. We have also computed N-L phase coexistence, but cannot 
resolve this (very concentrated, u ~ 3) regime accurately and do not present these results here. 

Can one observe l-L coexistence? This can only occur for concentrations above that necessary 
for equilibrium phase separation. One could conceivably prepare an equilibrium l-N mixture with 
the nematic phase in the log-rolling geometry. Upon applying shear, the system would then main- 
tain coexistence and move through the l-L two-phase region. However, the I phase is, itself, within 
the two phase region for l-N phase separation, so we expect the prepared coexisting l-L state to be 
metastable under shear. 

B. Common Strain Rate Phase Separation 



(a) (b) 
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u u 

FIG. 6. Common strain rate phase diagram in the (7— it) (a) and (a xy — u) (b) planes. The solid and 
dashed lines denote the phase boundaries for, respectively, the I and N phases. 

Common strain rate phase separation can be calculated exactly analogous to that for common 
stress phase separation. The resulting phase diagram for l-N coexistence, for L = 5, is shown in 
Fig. ^. The shear stress and composition are partitioned according to the lever rule in Fig. |5|b, with 

&xy = C&xyi + (1 - ()<Txyii, ( 2 -6) 

In this case tie lines connecting coexisting phases are parallel in the 7— u plane, and have a negative 
slope in the (a xy —u) plane because the I phase coexists with a denser and less viscous N phase. 
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There is an interesting crossover in the (a xy —u) plane. For dilute systems the stress in the N phase 
immediately outside the biphasic regime is less than the stress just before the system enters the 
biphasic region (Fig. |]b). Since the stress of the N branch is less than that of the I branch at the 
same strain rate and composition, we expect a decrease in the stress across the biphasic regime if 
composition effects are weak, e.g. near a critical point. For higher mean compositions the stress 
increases across the biphasic regime, because the width of the biphasic regime overcomes the 
shear thinning effect. 

Mean Constitutive Relations — The mean constitutive relation that could be measured in an ex- 
periment may be calculated from Fig. ^|a, and are shown in Fig. [7]. At higher concentrations the 
plateau has a positive slope while, coinciding with the crossover noted above, for lower concen- 
trations the plateau has a negative slope. A negative slope usually signifies a bulk instability, but 
here each band lies on a stable branch of its particular constitutive curve and the flow should be 
stable. Stable 'negative-slope' behavior has been seen in shear-thickening systems which phase 
separate at common stress [ |23p5| ], although in that case the mean constitutive curve was different, 
consisting of a backwards S curve and non-monotonic behavior only with multiple stresses for a 
given strain rate. 




0.0 0.2 0.4 0.6 0.8 0.0 0.1 0.2 0.3 0.4 0.0 0.2 0.4 0.6 

strain rate 7 * 10 4 7 * 10 4 7 * 10 

FIG. 7. Stress-strain-rate cr X y{i) curves for common strain rate coexistence. Same notation as Fig. ^ 

Based on analogies with equilibrium, we naively expect controlled strain rate and controlled 
stress experiments on concentrations such as those in Fig. [7]a-b to yield behavior similar to that 
for common stress phase separation, with nucleated or spinodal behavior depending on the applied 
strain rate, and the same caveat applying to decreasing the strain rate from agove. The situation for 
compositions with curves such as Fig. [7]a is qualitatively different. Here there is a range of stresses 
with three stable states: homogeneous I and N branches, and a banded intermediate branch. For 
controlled stress experiments, one possibility is that the I and N branches are favored in their 
respective domains of stresses. For example, in start-up experiments the system would remain on 
the I branch until a certain stress, at which point it would nucleate after some time. If the system 
nucleated onto the coexistence branch, increasing the stress further would return the system to the 
I branch. Since it nucleated from the I branch, it is more likely to jump directly to the N branch. 
Similar behavior is to be expected upon reducing the stress from the N phase. Another possibility 
is intrinsic hysteresis: that is, the system never jumps until reaching its limit of stability (from 
either the I or N side). The present theory cannot address this question. For controlled strain 
rate experiments, it could, in principle, be possible to maintain a stress on the two state region, 
although in practice this would also seem to be quite difficult, and would seem to be mechanically 
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unstable. In the case where stable composite curves with negative slopes were accessed, stress 
was the control variable 08[] . 



C. Common Stress or Common Strain Rate? 

What about the relative stability of phase separation at common stress or strain rate? While 
our one-dimensonal calculations cannot address this question, we have examined the two phase 
diagrams in the a xy — fx and j—fji planes, where // is the chemical potential QTOt] . This can be seen 
in Fig. |8]a-b where, for example, the I boundary for common strain rate phase separation (L ( ) lies 
in the N region for common stress phse separation, in the fi — a xy plane. This occurs because the 
stress of the I phase, at common strain rate, is larger than the stress of the N phase, due to the shear 
thinning nature of the transition. Conversely, the I phase at common stress lies within the I region 
of the common strain rate phase diagram. Analogy with equilibrium phase transitions suggests 
that, since the I phase of common strain rate phase separation thus lies on the "wrong" side of the 
phase boundary for common stress, given by the line in the [l — a xy plane, it would be unstable 
(or metastable) to phase separation at common stress. Conversely, the I phase at common stress 
is on the "correct" side of the coexistence line in the fi — 7 plane, and, again based on analogy 
with equilibrium, might be expected to be stable. Note that if the transition were shear thickening 
the situation would be reversed, and the arguments above would lead to common stress phase 
separation being unstable (or metastable) with respect to common strain rate phase separation. 

Boundary conditions may also play a role. In a Couette device the slight inhomogeneity of 
Couette flow induces an asymmetry between the inner and outer cylinders, exactly the symmetry 
of common stress phase separation (Fig. [I]). This should enhance the stability of common stress 
phase separation. Cone-and-plate rheometry induces a similar preference for the common stress 
geometry. 

(a) (b) (c) 
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FIG. 8. Phase diagrams in the (a) fi— 7 and (b) fi—o xy planes for I-N coexistence (the N state is stable for 
higher strain rate or stress, respectively). The thin vertical solid lines denote phase coexistence at common 
strain rate and stress in (a) and (b), respectively. The broken lines marked I 7 and iV 7 denote the coexisting 
states at common strain rate, in the Li—a xy plane (b); while the broken lines I a and N a , denote the coexisting 
states at common stress, in the li — 7 plane, (c) is the mean stress vs. strain rate curve. ABD denotes a path 
in the phase diagrams for common stress phase separation, while ABCD is a path which switches from 
common stress phase separation to common strain rate phase separation at B. 

An alternative possibility is presented in Fig. ||. If one argues that, in steady state, among 
the possible phases which are compatible with the interface solvability condition, the chemical 
potential reaches its minima so that no more diffusive material flux is possible. Based on such a 
criterion, upon increasing the strain rate for a given mean concentration the stable phase is that 
with the lowest chemical potential. The thick horizontal arrows in Fig. [8| denote the Li(cr xy ) and 
//(7) paths for the homogeneous high and low shear rate states, in the two phase diagrams. The I 
branch becomes unstable at A to phase separation at common stress, when the homogeneous path 
first crosses the phase boundary in the fx — a plane. For higher stresses the system follows the 
segment AB in Fig. [8|b, along the phase coexistence line at common stress, and follows the stress 
plateau AB in Fig. In the il — 7 plane the system phase separates, and the chemical potential 
as a function of mean strain rate follows the diagonal path AB in Fig. ||a (the dotted lines denote 
the strain rates of the coexisting phases). 

Upon increasing the strain rate further than point B, the chemical potential of the system can 
decrease by phase separating at a common strain rate. This reduces the chemical potential, at a 
given strain rate, from that of the segment BD to that of segment BC. Hence the system would 
take the path BC along the phase boundary in the /i — 7 plane, as far as point C, upon which the 
phase boundary crosses the homogeneous curve for the high strain rate phase of the given mean 
concentration. The path would be the diagonal path BC in the ii—o xy plane, and would correspond 
to the negative- sloped segment BC in the flow curve, Fig. |8]c. Finally the system follows the high 
strain rate branch, through CD. 

Upon increasing the controlled stress, the system would be expected to follow ABD. Upon 
decreasing the stress or the strain rate, DC-bottom jumping is expected. These scenarios follow 
from minimizing the chemical potential subject to the solvability constraint should phases coexist. 
Its correctness, of course, should be further examined by the full time evolution of the original 
dynamic equations. 

Experimental Studies — Mather et al. [|[ studied a liquid crystalline polymer melt (an aromatic 
polyester) and determined the lower limit of the l-N phase boundary in the 7 — T plane. The 
studies most relevant to the Doi model for rigid rod suspensions have been on wormlike micellar 
solutions near their isotropic-nematic coexistence region [|I],|39|], where common-stress banding 
was observed with a plateau stress that became steeper for concentrations closer to the equilib- 
rium l-N phase boundary, in qualitative agreement with our results. Common strain rate banding 
has not bee seen in these systems. Micelles are considerably more complicated than simple rigid 
rods, because they are not strictly rigid and their length (and hence coupling to flow) is a strong 
dynamic function of concentration. Experiments on micelles far from an apparent nematic tran- 
sition exhibit common stress shear banding with nearly flat coexistence plateaus P J^fniJ2Tj|r22| ], 



consistent with a concentration-independent JT^ ] instability (or transition). In kinetics studies the 
delay time before the transition to a banded (or high strain rate) flow in controlled stress start-up 
'quenches' diverged for a window of stresses slightly above the banding stress, whereas controlled 
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strain rate 'quenches' always decayed, eventually, onto a banded flow state. These interesting be- 
haviors cannot be explained by the topologies of the phase diagrams in Fig. |3|. Bonn et al. [27[] 
recently studied lamellar surfactant systems and observed slowly coarsening bands in the common 
strain rate geometry; and for controlled strain rate measurements they found transient constitutive 
curves analogous to Fig. [7]a or[7|b, consistent with common-strain rate phase separation. The true 
steady state behavior was not measured. 



IH. INTERFACE CONSTRUCTION 

Several microscopic and phenomenological models, as well as the apparent underlying flow 
curves for wormlike micelles, show an apparent degeneracy in the shear stress at which coexistence 
occurs (in the case of coexistence at a common stress). To resolve this degeneracy we have relied 
on the presence of inhomogeneous terms in the dynamical equations of motion, and determined the 
selected stress as that stress which allows a stable interfacial solution. In this section we explore 
this in more detail using a toy constitutive model. Similar arguments were given Spenley et al. 
[ ]33"| ] in a different languauge, and in a recent more rigorous study 

We consider planar flow with a velocity field v(r) = v(y)% with 7(3/) = dv/dy, and postulate 
the following constitutive relation for the shear stress: 

<r(i) = <rh(rr)-D(rr)8i<y. (3.1) 

The homogeneous flow curve 074(7) ls non-monotonic, as in Fig. ^|; and can be derived for a system 
with an underlying transition, as for the modified Doi model above, or from phenomenological 
models, such as the widely-used Johnson-Segalman (JS) model p^]. Gradient terms may come 
from the diffusion of the stress elements p3|]. The flow curve shown in Fig. ^|is for the JS model. 
In the (some what artificial) model where only the shear stress diffusion is considered, the steady 
flow condition for the JS model has the form of Eq. (P-ip, with D{^) oc 1/(1 + 7 2 ). 




1A IB 

applied strain rate 

FIG. 9. Non-monotonic homogeneous flow curve 0^(7) for the Johnson-Segalmann model. The thick 
curve denotes a banded flow between strain rates 74 and 7 S at stress gq. 

The steady state condition for planar flow is a uniform shear stress, 

a = 07,(7) - D(ff) d a 2 7 , (3.2) 
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with (T a constant. In an infinite system, an interfacial shear banding solution at a given stress <To 
satisfies Eq. <]3~2"1), with boundary conditions 



7(-ooj 
7(00) 



7a 
7b- 



(3.3a) 
(3.3b) 



9 v 7(±oo) = 



(3.4) 



Hence, given the second-order differential equation, the system is overdetermined. A solution is 
only possibly when these two conditions coincide, which may be obtained by varying the stress 
a - h is straightforward to integrate Eq. ( ^2| ) to show that a solution is possible when <t satisfies 
the following condition p5j] 



1A 



1A 



D(j) 



0. 



(3.5) 



Note that this is not an equal areas construction, unless D(y) is a constant D. 

Further insight may be obtained by casting the interface solution in terms of a dynamical 
system. Defining 



P = 7 

q = d y p= p, 



Eq. (3.2) becomes the following dynamical system, with y playing the role of time. 



p = q 

o"o - 07. (p) 



q 



D{p) 



(3.6a) 
(3.6b) 



(3.7a) 
(3.7b) 



For cr within the non-monatonic region of the flow curve the system has three fixed points 
P* = {Pa,Pb,Pc} on the axis q = 0, corresponding to the strain rates of the three homoge- 
neous flows. Linear stability analysis yields the stable and unstable manifolds of points A and B, 
with eigenvalues 



1 da h 
D(p) dp 



(3.8) 



p=p* 



and eigenvectors at angles 9 = arctan A with respect to the p-axis. Point C has imaginary eigen- 
values and is a cycle, while A and B are saddles with stable and unstable directions. 

An interfacial solution corresponds to an orbit connecting saddles A and B, and is denoted a 
saddle connection; it it also called a heteroclinic orbit, since it connects two different fixed points. 
This set of ordinary differential equations (ODE) does not generally have a saddle connection for 
an arbitrary a in the multi-valued region. It can be shown [ |35l ] that for models (with arbitrary 
numbers of dynamical variables) in planar shear flow with differential non-local terms that, apart 
from accidents, a saddle connection only exists at isolated points in the control parameter space. 



Here, the control parameters are cr and parameters which change the shape of cr^ Q35|J33Q , while 
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for the Doi model above, the control parameters are (for a given set of molecular parameters such 
as L) fi and a xy for common stress phase separation; and \i and 7 for common strain rate phase 
separation. 







j) 














FIG. 10. Orbits in the phase space for the dynamical system of Eqs. (p.7|). (Left) Away from the coexis- 
tence stress, a > (To, a saddle connection does not exist and A is connected to itself by a homoclinic orbit. 
[A similar situation exists for a < ctq, with the homoclinic orbit returning to B instead of A] (Right) At 
the coexistence stress a = a heteroclinic orbit connecting A and B (saddle connection) exists. The fixed 
points are at q = 0, and the solid lines show orbits beginning and ending on the stable fixed points A and 
B. 



Fig. [H] shows the evolution of "orbits" in the p — q phase space as the stress is tuned. For 
u = o"o a heteroclinic orbit exists, connecting A and B. This corresponds to an elementary shear 
band solution, in which one portion of the sample lies on the high strain rate branch B, another 
portion lies on the low strain rate branch A, and a single interface separates the two phases. For 
a ^ cr there is no heteroclinic orbit or saddle connection, and hence no stationary interface. 
Fig. [TO] (Left) shows a stress slightly greater than a , where a homoclinic orbit connecting state 
A to itself. Kramer []32|] pointed out in the context of reaction diffusion equations that such a 
homoclinic orbit corresponds to the critical droplet in a metastable phase of A material. Note that, 
although in real space it goes from A at y = — 00 to A at y = +00, the dominant spatial variation 
is in fact localized, with a size that vanishes when the stress reaches the maximum of the flow 
curve in Fig. ^ (at which the fixed points A and C annihilate). Slightly larger droplets are unstable 
and, when the full dynamics are returned to the problem, presumably flow to the high strain rate 
branch B, while smaller droplets are expected to decay back to A. By analogy with equilibrium 
behavior, for a > a we expect phase B to be the long time steady state, if fluctuations (i.e. noise, 
thermal or otherwise) were included. 



IV. CONCLUSION 

We have outlined the phenomenology of phase separation of rigid rod suspensions in shear 
flow, using the modified Doi model. Phase separation may occur with common stress or strain 



13 



rate, corresponding to the different coexistence geometry. We have calculated coexistence among 
three phases (paranematic, flow-aligning nematic, and log-rolling), while only two equilibrium 
phases exist. That is, the full rotational symmetry of an equilibrium nematic is broken by the 
biaxial shear flow, leaving two possible stable nematic orientations (the in-plane I and N states, and 
the out of plane L state). The shear thinning nature of the transition suggests that common stress 
phase separation is stable; while appealing to a minimization of the chemical potential, subjected 
to the interface solvability, predicts a curious crossover from common stress to common strain 
rate phase separation. We do not know which of these, or other, possibilities, are the physical 
ones. The composite stress strain curves depend on the coupling to composition [|23||, and can 
exhibit an apparent unstable constitutive relation, which would be mechanically unstable under 
controlled strain rate conditions. Although there have been few experiments on true lyotropic 
rigid rod systems in flow, wormlike micelles can have a flow-induced nematic phase at higher 
concentrations, and our results appear to qualitatively describe many aspects of these experiments. 
See, for example, the phase diagrams in Ref. [fHj]. 

We have also shown schematically how our construction for coexistence can be cast as an 
equivalent dynamical system, for which coexistence corresponds to a heteroclinic saddle connec- 
tion. In the most general case stress selection depends on the nature of the gradient terms in the 
dynamics [|3~5|], while in equilibrium systems the gradient terms can be exactly integrated to yield a 
condition independent of the gradient terms. The dynamical systems picture also yields an analogy 
with a critical droplet which may prove promising in understanding the non-equilibrium analogs 
of nucleation and growth. 
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